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Excitation spectra from angular momentum projection of Hartree-Fock states and the 

configuration-interaction shell-model 
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We make numerical comparison of spectra from angular-momentum projection on Hartree-Fock 
states with spectra from configuration-interaction nuclear shell-model calculations, all carried out 
in the same model spaces (in this case the sd, lower pf, and p-sd 5 / 2 shells) and using the same input 
Hamiltonians. We find, unsurprisingly, that the low-lying excitation spectra for rotational nuclides 
are well reproduced, but the spectra for vibrational nuclides, and more generally the complex specta 
for odd-A and odd-odd nuclides are less well reproduced in detail. 



O ■ I. INTRODUCTION 

(N 



In Hartree-Fock (HF) calculations, one approximates the ground state by a single Slater determinant (antisymme- 



U 

\Jl • terized product of single-particle wavefunctions) and minimizes the energy, reducing the many-body problem to an 
effective one-body problem [lj. The Hartree-Fock solution can break exact symmetries, such as rotational invariance, 
^,0 . and indeed is often more effective if it breaks exact symmetries. Restoring broken symmetries adds correlations that 
O^l ' further lower the energy and is a useful approach for many problems. 

In this paper we project states of good angular momentum from a Hartree-Fock state, and compare the resulting 
ground state energies and, especially, excitation spectra to exact results from equivalent full configuration-interaction 
diagonalization in a shell-model basis. Such comparisons help us to understand the accuracy and limitations of the 
projected Hartree-Fock approximation (PHF), as well as potentially providing a useful shortcut to nuclear structure. 
PHF is not new [2[ and can be found in a number of variants. Most of these are inspired by the Nilsson model 
[3, |4| and the Elliott SU(3) model [5j, and they can be classified by their use of (a) quasiparticle single-particle 
states (b) mixing of particle-hole states and (c) schematic or realistic interactions. For example, the projected shell 
model (PSM) uses quasiparticle single-particle states but no mixing of particle-hole states and schematic interactions 
Q; the MONSTER and VAMPIR codes use quasiparticle single-particle states, and mixing of particle- hole states, 
JTi | and realistic interactions [7j; and the projected configuration-interaction (PCI) eschews quasiparticles states but uses 
qv ■ realistic interactions and mixes higher-order particle-hole states Q. This list is not exhaustive. 

£SJ ' What is new in this paper is the comparison of exact numerical solutions against a relatively simple state, with a 

f — , focus on the quality of the excitation spectra. In this case the "exact" solution is from a full configuration-interaction 
_ C ' (CI) calculation, using semi-realistic interactions diagonalized in a truncated but nontrivial shell-model basis, against a 
single Hartree-Fock state with angular momentum projection. Hence we eschew quasiparticle states and particle-hole 
£f) • states, but we use a realistic interaction and we allow arbitrary deformation. 

This paper is similar in philosophy to previous work comparing the random phase approximation (RPA) against 
CI calculations in a CI basis [3. Ildj. testing how good is a specific approximation against a numerically exact result; 
in the former work the approximation was RPA, while here it is PHF. Another similar set of papers are those on PCI 
[8j], which compare results directly to CI diagonalization. An important difference with much previous work, however, 
is that in addition to studying the trends of the ground state energies, we consider in detail the excitation spectra. 

We find that rotational spectra, and spectra that are from simple particle-rotor coupling, are well reproduced, with 
a small rms error. Other excitation spectra, such as vibrational spectra or complex spectra from odd-A and odd-odd 
nuclides, have larger rms errors. This is unsurprising, but it is useful to verify in detail as we have done. 

In the next two sections we outline the numerical exact (CI) and approximate (HF and PHF) calculations, which 
we then compare for a variety of nuclides in Section V 

II. THE SHELL-MODEL BASIS AND CONFIGURATION-INTERACTION CALCULATIONS 

We begin by outlining the shell-model basis in which we work and the fundamentals of configuration-interaction 
calculations [Ill - [l3j . 

We work entirely in occupation or configuration space. This means that we start with some finite set of orthonormal 
single-particle states, {<p a (r)}. Each of these states have good angular momentum and parity and are labeled by orbital 
I, total j, and z-component j z ; implicitly they have good spin s, which here is 1/2 although other values are allowed 
by our code. We also have two species of particles, here protons and neutrons, each with fixed numbers. 



Using second quantization we have fermion creation and annihiliation operators (j>\,4> a that create or remove a 
particle from the states <p a (r). The uncoupled Hamiltonian is then 

ab abed 

The one-body and two-body matrix elements, T a b and V a b,cd respectively, are appropriate integrals over the Hamilto- 
nian and the single-particle states [111 ]. Because the Hamiltonian is an angular momentum scalar, we store the matrix 
elements in coupled form. Given the single-particle space, the matrix elements are generated externally to all our 
codes and read in from a file. Any radial form of the single-particle states and any form of the interaction, including 
non-local interactions, are allowed and make no pratical difference to our many-body codes. Because we work in 
second quantization, the two-body matrix elements are automatically antisymmeteric and we do not separate direct 
from exchange terms. 

For configuration-interaction (CI) calculations we use the BIGSTICK code [14] . BIGSTICK creates a many-body 
basis of Slater determinants {|a)}, constructed from the single-particle basis: 

i«>=n&io>, ( 2 ) 

a=l 

where N p is the number of particles (occupied states). We fix the total M — J z of the many-body basis states, 
(this is trivial as each single-particle state has good j z as well) and is thus called an M-scheme code; because the 
Hamiltonian is an angular momentum scalar, the final eigenstates thus will automatically have good total angular 
momentum J (and isospin T, although we have the capability to break isospin; T z = (Z — N)/2 is fixed). Within the 
fixed single-particle space and fixed J Zl T z we allow all possible configurations. 

Given the basis, BIGSTICK then computes the many-body Hamiltonian matrix elements (a\H\/3) and finds the 
low- lying eigenstates, including the ground state, using the Lanczos algorithm [13, [l5[- With the single-particle space 
and Hamiltonian matrix elements fixed, the resulting eigenenergies for the full-space CI calculation are numerically 
exact. BIGSTICK can handle CI M-scheme model spaces on a desktop computer up to dimension roughly 2-400 
million. 

III. HARTREE-FOCK APPROXIMATION 

The Hartree-Fock approximation (HF) is a variational method using a single Slater determinant [1]. One applies a 
unitary transformation among the single-particle states: 

N e 

Ct = ^U ai (j> a , (3) 

o=l 

where N s is the number of single-particle basis states (hence the number of particles N p < N s ) and creates the trial 
Slater determinant 



Then one finds a Slater determinant that minimizes the energy, that is, that minimizes 

S = <*rlfl*f> (5) 

Because our input matrix elements are antisymmeterized, we fully include both direct and exchange terms. Working 
in occupation space the exchange term causes us no difficulty (or, to put it another way, any difficulty is off-shored 
into calculation of the antisymmeterized integrals) . 

Many Hartree-Fock calculations enforce good angular momentum, for closed-shell systems or closed-shell plus or 
minus one particle. Our HF code (the SHERPA code Q) allows the Slater determinant to break rotational invariance, 
even for closed-shell systems; such Slater determinants are 'deformed.' Experience in nuclear physics suggest this to 
be a fruitful path. Our only constraint is that we assume the transformation ([3]) to be real. 



In order to minimize, we use the standard Hartree-Fock equations, which consist of iteratively solving 

hili = e»Ui, (6) 

where the Hartree-Fock effective one-body Hamiltonian is 

h a b = T a b + 2_^ Vac,bdPdc (7) 



cd 



and 

Pdc = Yl U ^d (8) 



JV 



We have separate Slater determinants for proton and neutrons; the generalization is straight forward. More details 
are found in Appendix [X] 

IV. PROJECTION OF ANGULAR MOMENTUM AND PARITY 

In order to project out angular momentum, we introduce the standard projection operator [l|, |4| 

P J MM , = / <m-D^,(Sl)R(Sl), (9) 



where O = a, /3, 7 and dfl = da sin (3 d (3 'dj, 'D M \ I , is the Wigner D-matrix [l(j, and R(Q) is the rotation operator. The 
rotation operator acts on a Slater determinant U whose columns are single-particle states with good j, m; therefore 
the matrix elements of R are given by the Wigner Z?-matrix: 

(j'm'\R\jm) = Sj.jiy^JaM)- (10) 

It is useful to note that the j, to, to' for the rotational matrix R are those of the single-particle space, while the 
J, M, M' for the projection operator (j9]) are those of the many-body space. 
We now introduce the Hamiltonian and "norm" matrix elements, 

h J M>M , = {*\HP J MtM ,\*), n J MM , = mP J MtM ^) (11) 

These are both Hermitian. 

Then we solve the generalized eigenvalue equation 

Yl h M,M>9M> = EJ Y. n M,M>9 J M>- (12) 

M' M' 

where we allow qm to be complex. Although the deformed Hartree-Fock state will have an orientation, the final result 
will be independent of orientation; we confirmed this by arbitrarily rotating our HF state. 

PHF will only generate a limited number of states. For a given J, the maximum number is (2 J + 1), although 
fewer can be found if the norm matrix n has zero, or very small, eigenvalues. This can happen, for example, if the 
HF state has axial symmetry. Therefore there will be states missing from the CI spectrum, such as low-lying excited 
J = states. 

V. RESULTS 

In this section we apply our calculations to a number of nuclides in the sd and pf shells as well as a limited p-sd 
space where we allow for parity-projection as well. We ask a number of questions: 

• Because PHF is a variational theory, the PHF ground state energy will be above the exact ground state energy. We 
ask: what is the systematics of that displacement? For example, if the displacement were nearly constant, we could 
use PHF to reliably estimate the true ground state energy. 

• In a similar fashion we ask, how well does PHF yield the excitation spectrum? 

• For odd-A and odd-odd nuclei, how often does PHF yield the correct J w assignment for the ground state? 



• For multi-shell calculations, how well does parity+angular-momentum projection yield the splitting of positive and 
negative parity spectra? 

Many previous studies have focused on the ground state energy. By contrast, we focus on the excitation spectrum, 
as well as systematics of the ground state energy. 

We work with several model spaces and interactions. All model spaces assume some inert core and valence particles; 
the single-particle states can be thought of as spherical harmonic oscillator wavef unctions, but that has no impact 
on our calculations. The interactions we use are all 'realistic,' based on effective interactions derived from scattering 
data but with individual matrix elements adjusted to fit many-body data (see [111 Il7j for details of methodology). 
The three model spaces we work in and their interactions are: 

• sd, or the ls -i /s- Oria/s-Ody? valence space, assuming an inert 16 core; the interaction is the universal sd-interaction 
'B,' orUSDB [l3; 

• pf, or the IP1/2-IP3/2-O/5/2-O/7/2 valence space, assuming an inert 40 Ca core; the interaction is the monopole- 



modified Kuo-Brown G-matrix interaction version 3G, or KB3G 18[; 
• and finally the psd, or 0p3/2-0pi/2-0c?5/2- lsi/2 valence space, assuming an inert 4 He core; the interation is a hybrid 
of Cohen-Kurath (CK) matrix elements in the Op shelipjj, the older universal sd interaction of Wildenthal [20] in 
the 0d 5 / 2 -lsi/2 space, and the Millener-Kurath (MK) p-sd cross-shell matrix elements [211] . We leave out the 0d 3 / 2 
orbit is to make full shell-model calculations tractable. Within the p and sd spaces we use the original spacing of 
the single-particle energies for the CK and Wildenthal interactions, respectively, but then shift the sd single-particle 
energies up or down relative to the p-shell single particle energies to we get the first 3 _ state at approximately 6.1 
MeV above the ground state. The rest of the spectrum, in particular the first excited + state, is not very good, but 
the idea is to have a non-trivial model, not exact reproduction of the spectrum. This model space and interaction 
allows us to consider model nuclei with both parities and to investigate parity-mixing in the HF state. 

To illustrate our results, we begin with 24 Mg in the sd shell. Fig. Q] compares the exact CI spectrum on the right 
against the PHF spectrum on the left. The PHF g.s. is 1.78 MeV above the CI g.s., while on average each states 
in the PHF spectrum is shifted by 1.68 MeV. For comparison, the correlation energy, here defined as the difference 
between the unprojected HF energy and the CI g.s. energy, is 6.62 MeV, so that restoring good J by projection 
accounts for 73% of the correlation energy. 

Fig. [2] is the same as Fig. [1] but with the PHF spectrum shifted down so the g.s. energies coincide. The agreement 
between the excitation spectra is very good, with an rms error of 0.10 MeV. 

Other nuclides shall qualitatively similar results, albeit with varying degrees of accuracy. In general the PHF 
excitation spectra of even-even nuclides was of much higher accuracy than for odd-odd and odd- A nuclides, which is 
not too surprising as one might expect a deformed HF state to approximate a rotational band. 

For example, consider the odd-odd nuclide 30 A1 in the sd shell, as shown in Fig. [31 The PHF g.s. is 6.40 MeV 
above the CI g.s., and PHF only accounts for 25% of the correlation energy. When the PHF g.s. is shifted down to 
the match the CI g.s., Fig. [H the rms error in the excitation spectrum is 0.66 MeV and many of the states are out of 
order-in particular the first excited 1 + state. 

This continues in the pf shell. In Fig. [5] we show 52 Ti, which has a strong vibrational rather than rotational 
spectrum. Here the PHF g.s. energy is 3.16 MeV above the CI g.s., account for only 42% of the correlation energy, 
and even when the PHF spectrum is shifted down in Fig. [5J it is clearly compressed relative to the exact CI, with an 
rms error of 0.53 MeV. 

By contrast, the spectrum of 49 Cr, which can be thought of as a particle-rotor coupled nucleus, is good. As shown 
in Fig. [71 the PHF g.s. is 2.63 MeV above the CI g.s., and PHF accounts for only 35% of the correlation energy; but 
when the PHF spectrum is shifted down, as in Fig. [8] the rms error is only 0.01 MeV. 

We also worked in a cross-shell system, the p-sd 5 / 2 space. Keep in mind we artificially reduced the separation 
between the p and sd 5 / 2 shells in order to provoke a parity-mixed HF state. 

The vibrational spectrum the positive parity states of 18 C is only modestly reproduced by PHF, as illustrated in 
Fig. [91 andfTOl The PHF g.s. is 3.84 MeV above the CI g.s., accounting for only 44% of the correlation energy, and 
when the PHF spectrum is shifted to match the g.s. energies, the rms error is 0.22 MeV, with the 4 + state in the 
wrong place. 

Somewhat better are the positive parity states of 22 F, in Figs. [TT1 and [T2l The PHF state is 1.41 MeV above the 
CI g.s., but that only accounts for barely 15% of the correlation energy. On the other hand, when the g.s. energies 
are matched, the error in the excitation is only 0.08 MeV, although some of the states are in the wrong order. 

Tables I and II summarize our results. In particular, Table I summarizes for specific cases the standard deviation 
(rms error) of the PHF excitation spectra relative to the CI excitation spectra, once the ground state energies are 
matched. Table II summarizes all our calculations; the rightmost column gives how often the correct g.s. ,P arises 
out of the PHF. Even-even nuclides almost always get the correct J = ground state (but for our mixed-shell 
calculations, occasionally the wrong parity-but keep in mind we artificially reduced the splitting between the shells 
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FIG. 1: Comparison of exact CI and PHF spectra of Mg in the sd model space with the USDB interaction. 
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FIG. 2: Comparison of exact CI and PHF spectra of 24 Mg in the sd model space with the USDB interaction, but with the 
PHF spectrum shifted downwards so the g.s. energies coincide. 

to force a mixed-parity HF state), but odd- A got the correct g.s. assignment only a little over half the time, and 
odd-odd nuclei less than half the time. 



VI. CONCLUSIONS 



We have carried out detailed angular momentum projected Hartree-Fock calculations in a shell-model basis with 
semi-realistic shell-model interactions, and compared the spectra to exact configuration-interaction diagonalization 
calculations, with a particular focus on the quality of the excitation spectra. We found, unsurprisingly, that rotational 
nuclei are best reproduced by PHF. 

The U.S. Department of Energy supported this investigation through grant DE-FG02-96ER40985. 
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FIG. 3: Comparison of exact CI and PHF spectra of Al in the sd model space with the USDB interaction. 
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FIG. 4: Comparison of exact CI and PHF spectra of Al in the sd model space with the USDB interaction, but with the PHF 
spectrum shifted downwards so the g.s. energies coincide. 



TABLE I: Summary of data for the nuclei in the paper. The standard deviation is for the 
shifted PHF spectra relative to the exact CI spectra. 



Shell 


Nucleus 


Valence 
Z N 


Std. dev. 

(MeV) 


% PHF 
states 


sd 


24 Mg 
30 Al 


4 4 

5 9 


0.077 
1.004 


7 
12 


Pf 


52 Ti 
49 Cr 


2 10 
4 5 


0.714 
0.074 


8 
4 


P-sd 5/2 


18 C 
22p 


4 10 
7 11 


0.810 
0.194 


7 
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FIG. 5: Comparison of exact CI and PHF spectra of 5 Ti in the pf model space with the KB3G interaction. 
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FIG. 6: Comparison of exact CI and PHF spectra of 52 Ti in the pf model space with the KB3G interaction, but with the PHF 
spectrum shifted downwards so the g.s. energies coincide. 



Appendix A: Some computational details for Hartree-Fock 



A single general Slater determinant we write as a rectangular, N s x N p matrix \&. Operationally, the vectors Ui, of 
length N s , from solving the Hartree-Fock equation ^ form the columns of the transformation of the single-particle 
basis ([3]). We take the N p vectors with lowest e\ to form the columns of *& . In other words, *f? is simply selecting 
columns of U, which are computed in Eq. (J6j> ; this is the unitary transformation that diagonalizes h . We have 
separate proton and neutron Slater determinants, ^^ and M/,,. 

From the Slater determinant we construct the one-body density matrix, 



N„ 



Pab = ( * 



ftj b *)=£* ai *k. 



i=l 



We then construct the Hartree-Fock effective one-body Hamiltonian, 



h a b — H ab + 2_^ V ac ,bdPcd- 
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FIG. 7: Comparison of exact CI and PHF spectra of Cr in the pf model space with the KB3G interaction. 
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FIG. 8: Comparison of exact CI and PHF spectra of 49 Cr in the pf model space with the KB3G interaction, but with the PHF 
spectrum shifted downwards so the g.s. energies coincide. 



Because we have two species, we have to take a little care. If our interaction is an isospin scalar, then 



v pp = v nn = V T=1 

abed abed v abed 



and 



(A3) 



vT; = - (v T h =? + v T t r}) 

abed n V abed ' v abed ) 

Here we have suppressed the obvious decoupling of angular momentum [111 . Il6| ] . 
Then the proton Hartree-Fock Hamiltonian is 

h 7 " — H 1 -4- V^ V pp n^ 4- V pn n v 
n ab — n ab T 2_^ V ac,bd Pcd + V ac,bdP cd 

cd 



(A4) 



(A5) 



and similarly for the neutron Hartree-Fock Hamiltonian. 
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FIG. 9: Comparison of exact CI and PHF spectra of C in the p-sd 5 /2 model space. 
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FIG. 10: Comparison of exact CI and PHF spectra of C in the p-sd 5 / 2 model space, but with the PHF spectrum shifted 
downwards so the g.s. energies coincide. 



Appendix B: Some computational details for projected Hartree-Fock 

We now have to compute (^>\R(^l)\ip) and fi>\HR(£l)\t[j), where R is the rotation operator. In our case 

R\i&) -> R* = * 



(Bl) 



where ^ and >& are N s x N p matrices representing our original and rotated Slater determinants, respectively, and R 
is the square matrix the carries out the rotation. In our single-particle basis the matrix elements are straightforward: 



Rab = (jam a \R(a,P 1 j)\j b m m ) = 5 jajh T>% (a,/3,j), 



(B2) 



where V 3 m , m is the Wigner D-matrix [H 



Compute the matrix elements between two oblique Slater determinants is relatively straightforward [22j. First, 

(#|£|tf} = (*|#) =det* f ^. (B3) 
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FIG. 11: Comparison of exact CI and PHF spectra of F in the p-sds/2 model space. 
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FIG. 12: Comparison of exact CI and PHF spectra of F in the p-sd 5 / 2 model space, but with the PHF spectrum shifted 
downwards so the g.s. energies coincide. 



Calculation of the Hamiltonian (and other) expectation value requires the density matrix, 



Pab = 






* 



(V*) *t 



ba 



Then 



(*|iJ|#) = J2 H lbP'ab + Yl V *b,cd(PacPbd - PadP'bc 



(B4) 



(B5) 



'ibcd 



Accounting for two species is straightforward. Note that even if our original Slater determinants were real, by rotating 
over all angle they can become complex. 

We only project out good angular momentum. We could in principle project out exact isospin; to do this we would 
have to use a single Slater determinant containing both protons and neutrons and rotate in isospin space. We leave 
such a modification for future work. 



11 



TABLE II: Frequency of ground state matching. The "No. nuclei" column represents 
the total number of nuclei that was studied for that shell that satisfies the Z-N 
relationship. 







No. 


Correct 


Incorrect 


Frequency 


Z - N 


Shell 


nuclei 


G.S. J 


G.S. 


J 


correct 




s-d 


15 


15 







1.00 


Even-Even 


P-f 


10 


10 







1.00 




p-sd 


17 


12 


5 




0.71 




s-d 


15 


8 


7 




0.53 


Odd-Odd 


P-f 


6 


2 


4 




0.33 




p-sd 


12 


4 


8 




0.33 




s-d 


25 


17 


8 




0.68 


Odd-A 


P-f 


11 


6 


5 




0.55 




p-sd 


25 


13 


12 




0.52 



Appendix C: Linear algebra of projected states 

In order to calculate the weights of the angular momentum-projected states in the Hartree-Fock state, we have to 
investigate in some detail the linear algebra of solving Eq (|12p . 
Specifically, we want to expand the Hartree-Fock state 



|#F)=$>a|A>. 

A 

Now each state A has definite J, but can be broken up into M-components. Thus we expand 



X,M 



and we ultimately want 



«A 



M 



a\.M\ 



As discussed above, we want to solve the generalized eigenvalue problem 

Ugx = E x Afgx 
Now because 

Y d JKMP J K<M = l 



(CI) 



(C2) 



(C3) 



(C4) 



(C5) 



we have txM — 1. 

To solve the generalized eigenvalue problem we decompose the norm matrix, J\f — SS\ for example by Cholesky 
or by diagonalization. In general there will be a null space which we project out. 

By creating h = S~ 1 HS~ 1 ' ! (and also projecting out the null space, to eliminate division by zero), we then solve 
the ordinary eigenvalue problem 



h^A = E x v\. 



(C6) 



Of course, \v\\ = 1. Furthermore, because of this, the vectors g\ = mathbfS ^v\ are orthonormal with respect to 



the norm matrix, that is 



Another way to see this is to note that 



51-AAsp 



3a m - 



\X,JM) = J2Pm,k\HF) 9k (X) 



(C7) 
(C8) 



K 
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and because P 2 — P we automatically get (|C7|) . 

As stated above, our ultimate goal is to compute \a,\\ 2 - Towards this end, we note that 



so that 



a x ,M = (A, JM\HF) = J2^m,k9kW (C9) 

K 



\ax\ 2 = J2 I«a,m| 2 = E te(A)(AA 2 )L x , 9 ^(A) = .9|AA 2 gA. (CIO) 



M if, if' 

We can rewrite this as 

\a x \ 2 =v{S^Sv x (Cll) 

which is what we want in convenient terms. To confirm this is correct, we compute 



A A V A / 



(C12) 



but using the completeness relation, the sum over A yields just 1, and we have 

tr55 t = trAT=l (C13) 

as we'd expect. 
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